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ABSTRACT 



Completion of the installation at the Naval Postgraduate 
School of a cold-flow test facility for the turbine of the 
Space Shuttle Main Engine High Pressure Fuel Turbopump is 
reported. The article to be tested is the first stage of 
the 'Alternate Turbopump Development' model designed and 
manufactured by Pratt & Whitney. The purpose of the 
facility is to enable the development of non-intrusive flow 
measurements and comparison of those measurements with 
numerical sinulations. Flow field characteristics of the 
turbine stator were predicted using a three-dimensional 
viscous flow code. A sensitivity study was conducted to 
determine the effect of inlet profile to flow field 
solution. Recommendations are made for future use of the 
test facility and validation of the numerical simulation 
scheme . 
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I . 



INTRODUCTION 



A. PURPOSE 

This thesis describes the modification, completion and 
initial test of a test facility designed to enable the 
development of non- intrusive measurement techniques within 
the rotating blade row of a transonic turbine. It also 
details the computer modeling of a turbine stator flow field 
for validation of an advanced viscous flow code. 

B . OVERVIEW 

The performance of a turbomachine depends heavily on the 
flow fields generated within the machine during it's 
operation. Measurement, prediction and understanding of 
conditions in the flow field are essential to optimized 
design of turbomachines. More efficiently cooled blading, 
enabled by advanced computer modeling and design, can result 
in significant improvements in service life, reliability and 
lifecycle cost . 

Installation at the Naval Postgraduate School 
Turbopropulsion Laboratory (TPL) of the cold- flow test 
facility reported by Studevan [Ref. 1] was completed as part 
of the present work. The facility was designed to enable 
tests of a turbine intended to power the High-Pressure Fuel 
Turbopump (H TJ FTP) (Fig. 1) of the Space Shuttle Main Engine 
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Figure 1 . Space Shuttle Main Engine 
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(SSME) . The turbine was designed and manufactured by Pratt 
& Whitney as the 'Alternate Turbopump Development' (ATD) 
model . 

The purpose of the new Turbine Test Rig (TTR) 
configuration (Fig. 2) at TPL is to facilitate the 
advancement of non- intrusive Laser -Doppler Velocimetry (LDV) 
techniques required to validate viscous flow simulation 
codes. The long term goal of this research is to use LDV to 
measure the flow field immediately adjacent to the endwall 
of the turbine in the tip clearance region. An intermediate 
goal is to determine the performance of the single, and 
ultimately the two- stage turbine operating continuously in 
the TPL facility. Gaddis, Hudson and Johnson [Ref. 2] 
documented the exit conditions and overall performance of 
the ATD, while Hudson and others [Ref. 3] investigated the 
original turbine model. Both of these tests were conducted 
in a short drration blow-down facility at NASA Marshall 
Space Flight Center . 

1. Experiment 

Initial measurements of the ATD will determine the 
power output of the single stage turbine. The 
instrumentation documented herein will be primarily that 
concerned with the safe operation and control of the 
facility. Provisions have been made for initial flow field 
testing, which will be made with classical intrusive 
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techniques. The initial tests were conducted to determine 
the integrity of the turbine and its subsystems. 

2 . Numerical Simulation 

Numerical simulations of the flow field in the first 
stage stator of the SSME HPFTP turbine were accomplished 
using two codes developed by Roderick V. Chima called 
Turbomachinery C-Grid (TCGRID) [Ref. 4], and Rotor Viscous 
Code 3-D (RVC3D) [Ref. 5] . TCGRID was used for the 
generation of the 3D volume grid around a single blade 
element and RVC3D performed the numerical solution of the 
f lowf ield . 
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II. FACILITIES 



A. EXISTING TEST CONFIGURATION 

The existing configuration of the TTR, as reported by 
Studevan [Ref. 1], consisted of: 

• Inlet Ducting 

• ATD 

• Shaft and Bearing Housing 

• Bearing Housing Support 

• Dynamometer 

• Bearing Lubrication System' 

• Controls and Instrumentation. 

All of the above listed components are contained within the 
TTR cell with the exception of the Controls and 
Instrumentation which are housed in a control room adjacent 
to the cell. The interior of the TTR cell may be seem from 
the control room through an explosionproof window. This 
allows the test rig to be operated and observed from a safe 
environment. The test arrangement is shown in Figure 2. 

During operation, high pressure air flows from the inlet 
ducting through the inlet strut housing and into the ATD. 

In the ATD, air passes through the annulus of the first 
stage stator where it is accelerated and turned off axis. 

The air then impinges on the first stage rotor blades which 
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move in reaction to the pressure, turning the rotor disk and 
converting the energy in the air to rotational energy. The 
air then exhausts into the TTR cell. As the rotor disk 
turns, it spins the main shaft in the bearing housing. The 
main shaft, via the quill shaft, is coupled to the 
dynamometer where the rotational energy is absorbed. 

Shaft speed of the turbine may be sensed and controlled 
by the dynamometer system. Configuration of the speed 
control system is unchanged from that reported by Kane [Ref. 
6] . Speed may be automatically or manually controlled by 
adjusting the water flow rate through the dynamometer . In 
the automatic mode, a magnetic speed- sensor signal is 
compared to a set point, and a controller drives the 
difference to zero. The dynamometer is instrumented to 
provide mechinical power output via torque measurements. 

The power output may be calculated in two other ways; 
namely, using the mass flow rate and temperature rise of the 
water through the dynamometer, or using the mass flow rate 
and stagnation temperature drop of the air through the 
turbine . 

B. CONFIGURYTION CHANGES 

1 . Bearings 

The original design of the shaft and bearing housing 
called for four high-precision ball-bearing units. After 
coupling the shaft and bearing housing with the ATD, the 
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assembly was to be dynamically balanced. Balancing proved 
to be impossible and the assembly was taken apart for 
inspection. Upon disassembly, one of the bearing units was 
found to be damaged. In order to simplify the system, the 
outer bearing spacer was redesigned so that only two ball 
bearing units were necessary (Fig. 3) . Appendix A contains 
the engineering drawing. 

Duriig reassembly of the shaft and bearing housing, 
the shaft was cooled in a refrigerator to 62 'F and the 
bearings were slowly heated in an oil bath to 163 ’F. It was 
noted that condensation formed on the shaft when it was 
removed from the refrigerator 'and it was wiped down prior to 
inserting the bearing. The shaft was placed on end on a 
workbench. When the bearing reached the proper temperature, 
it was removed from the oil bath with tongs and placed on a 
metal plate. From the plate, it was immediately picked up 
in gloved hands and carefully but quickly inserted over the 
shaft. When the bearing was inserted over the shaft in an 
attitude level with the workbench it did not bind and slid 
properly into position. The procedure was practiced several 
times with the damaged bearing before an attempt was made 
with the actual hardware. To insert the shaft and bearing 
assembly into the bearing sleeve, the sleeve was preheated 
on a hot plate to 100 ‘F. The shaft and bearings were then 
pressed into the bearing sleeve with a milling machine used 
as a press. 
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Figure 3. Bearing and Shaft Showing Oil Mist Path 
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Oil mist sump 



2. Oil-Mist Cooling System 

The Portable Turbine Lubrication Unit described by 
Studevan [Ret. 1] for cooling the high-precision ball 
bearings was replaced by an oil -mist lubrication system. 

The system was powered by shop air which passed through a 5 
micron filter to remove particulate matter before entering 
the Norgren Micro-Fog Unit, where it forced a small stream 
of atomized oil through tubing to a nozzle (Fig. 4) . The 
nozzle was inserted into a hole atop the downstream end of 
the bearing sleeve (Fig. 3) . The oil/air stream flowed 
through openings in the outer bearing spacer and impinged on 
the ball bearings. It then passed through the bearings and 
out two ports on the lower upstream and downstream ends of 
the bearing sleeve. This was an open system. The used oil 
vented to waste buckets and was not recirculated. 

3. Temperature Monitoring System 

Because inadequate cooling could lead to bearing 
degradation, temperature - sensing thermocouples were placed 
in holes through the bearing housing at each bearing outer 
race. These were wired to digital output meters installed 
in the control panel in the TTR control room and were used 
to constantly monitor the bearing temperature and insure 
that the cooling system was operating effectively. The 
redesigned bearing housing assembly was rebalanced 
successfully at approximately 2500 rpm , with the oil -mist 
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Figure 4. Oil Mist Cooling System 
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lubrication oystem operating. The temperature-monitoring 
system was installed, and the assembly was put into service. 

4. Vibration Monitoring System 

In order to ensure the safety of operating the TTR, 
a vibration - detection system was installed to detect 
possible out - of -balance conditions and bearing degradation. 
The vibration - sensing system consisted of four 
accelerometers and a signal analyzer . All four 
accelerometers were attached to the bear ing - housing support. 
One was upstream at the top of the housing and the second 
was upstream but displaced 90' clockwise looking upstream. 
The remaining two accelerometers were installed on the two 
bearing -housing support struts. The first two were wired 
directly into a Scientific -Atlanta SD380 Signal Analyzer in 
the control coom. The two accelerometers installed on the 
support struts were wired into the control panel, and to an 
analog Ballantine Laboratories RMS indicator . From the 
analog indicator, they could be connected (if needed) to the 
signal analyzer for digital processing and hard copy output. 

Prior to powering the turbine with air flow, it was 
analyzed to determine natural frequencies so that resonance 
at those frequencies could be avoided. The accelerometers 
were installed and wired to the signal analyzer before any 
of the inlet ducting was attached. The stationary turbine 
and bearing housing assembly was struck with a modally tuned 
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hammer and the response of the accelerometers was recorded 
by the signal analyzer . An example from the top 
accelerometer may be seen in Figure 5. The top trace shows 
the real time response and the bottom shows the frequency 
spectrum with the first predominant mode at 67.5 Hz. The 
assembly was then turned by electric motor and a belt drive 
at varying speeds up to 960 rpm and the largest amplitude 
responses were found and graphed, as shown in Figure 6. 

Trace A is from the top accelerometer and shows the 16 Hz 
speed signal along with its higher harmonics. Trace B is 
from the accelerometer at 90', and shows similar 
characterist lCs . 

5. Control System 

The control system in the TTR control room included 
a Fisher - Por ter electronic set-point controller for the 
outlet valve of the hydraulic dynamometer . This was used as 
a speed regulator with manual or automatic modes and was 
unchanged from the configuration described by Studevan [Ref. 
1] and Kane [Ref. 6] . In the automatic mode, the controller 
received a conditioned feedback signal from a magnetic speed 
sensor in the dynamometer. The controller compared the 
speed signal to a set-point and produced a corrective 
signal. The corrective signal, received by the outlet 
valve, controlled the water level in the dynamometer. The 
speed varied in response to the resistance produced by the 
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Figure 6. Motor Actuated Rotations] Vibration Response 
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water and the energy was dissipated into the water as heat. 
The manual inlet control valve was also located in the TTR 
control room. Manual operation of the speed control system 
was used during all the tests conducted in the present 
study . 

6 . Probe Placement 

A probe access hole was cut through the outer casing 
of the turbine upstream of the first-stage stator. The hole 
location corresponded with the inlet of the grid designed to 
be used in the computer simulation of the flow field. A 
probe holder was designed and manufactured to hold the probe 
in place on the test rig. It was designed to be used in any 
position around the circumference of the turbine. However, 
a new access hole would need to be cut in order to use it in 
another position. An engineering drawing of the probe 
holder may be found in Appendix A. 

CombinatLon probes (numbers 2 and 3) were placed in the 
inlet duct through holes cut in the piping. They were 
positioned at the same axial station to measure stagnation 
temperature and stagnation pressure into the turbine. 
Combination probe number 4 was clamped to the table 
immediately downstream of the turbine to measure exit 
stagnation temperature and pressure. 
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Ill . EXPERIMENT 



A . EXPERIMENTAL PROCEDURE 

Before any measurements were attempted in the TTR, the 
shroud inserts between the blade tips and the turbine outer 
casing were removed. With the shroud inserts removed, there 
was no danger of blade tips rubbing due to misalignment. 
Initial system validation runs documented herein were 
conducted with the shroud inserts removed. The TTR was 
operated twice to gain experience, test safe operating 
methods and shake-down' the facility. During the first 
shake-down run, speed control and vibration response were 
tested under various air and water flow settings. It was 
noted that if water flow was sufficiently reduced by closing 
the outlet valve to the dynamometer, the speed of the 
turbine climbed rapidly. On inspection of the TTR after the 
run, the dynamometer was found to be very hot. During 
inspection of the water - circulation system, water was found 
to have overflowed onto the floor of the pumping station due 
to surging developed by changing demand to the dynamometer . 
It was decided that the water flow through the dynamometer 
should be set at the start of the run to keep the 
dynamometer cool and eliminate surging of the water 
circulation system. 
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For the second shake -down run, a specific procedure was 
adopted to set the water and air flow at the start of the 
run. In order to keep sufficient and consistent water flow 
through the dynamometer, the inlet valve to the dynamometer 
was set to 80% open and the outlet valve was set to 50% 
open. For mass flow to the TTR to be properly calculated, 
all the air flowing through the metering orifice had to go 
to the turbine (Fig. 7) . To ensure that all air flow went 
to the turbine, the number - two dump valve had to be 
completely closed. After the TTR inlet valve was set to 20% 
open, the number - two dump was closed, then the number -one 
dump valve was set to 40% open: Further manipulation of the 

mass flow to change the speed of the turbine involved 
opening the TTR inlet valve to increase speed or opening the 
number -one dump to decrease speed. It was necessary to have 
some air venting through dump valve number one in order to 
keep the total flow rate sufficiently high or the Allis- 
Chalmers compressor (Fig. 7) could surge due to excessive 
back pressure. 

During the experimental run, the procedure described 
above for setting water and air flow was followed. Pauses 
were made during the acceleration and system settings were 
noted. At each pause, the speed was allowed to stabilize 
and measurements were made of temperatures and pressures for 
the power output and mass flow calculations. Measurements 
relevant to power calculations are shown in Table I . 
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TABLE I . ABBREVIATED TTR DATA 



Speed 


ti 


K 


Pi 


TT4 


TT3 


(RPM) 


CF) 


( in . H z O) 


(in. Hg) 


CF) 


CF) 


880 


81 


0 


37 .35 


77 


83 


1420 


88 


19 . 0*9 


38 . 57 


83 


88 


1560 


92 


19 . 000 


38.98 


86 


92 


2100 


94 


16 .625 


41.42 


'll 


33 


3700 


97 


24 . 500 


53.23 


84 


97 


5950 


105 


33 .750 


49 . 56 


36 


105 


7900 


111 


17.375 


55.06 


71 


111 



B. MASS FLOW AND POWER CALCULATIONS 

The turbine air mass-flow rate was metered by an orifice 
plate (Fig. 7) and was calculated using the method described 
by Vavra [Ref. 16] and summarized in Appendix B. The power 
output of the turbine was calculated using the stagnation 
temperature drop of the air across the turbine, and the mass 
flow rate. Table II shows calculated values of mass -flow 
rate, specific power out of the turbine, and efficiency. 
Details of these calculations may be found in Appendix C and 
a complete table of measurements is contained in Appendix D. 
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TABLE II. TURBINE PERFORMANCE (SHROUDS REMOVED) 



1 

Speed 

(RPM) 


Mass Flow 
Rate 
(lbm/s) 


Specific 

Power 

(Hp/lbm) 


= ^ i 

Efficiency 

(%) 










880 


- 


- 


14 .96 


1420 


7 .51 


1 .35 


7 .47 


1560 


7 .63 


2 . 03 


9.98 


2100 


7 . 37 


2 . 36 


6 .48 


3700 


10 . 08 


4 .40 


3 . 86 


5950 


11 . 27 


9 .81 


3 .49 


7900 


8.57 


13 . 53 


3 .23 
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IV. NUMERICAL SIMULATION 



A. GRID GENERATION 

Numerical simulation of a flow field is highly dependent 
on the grid used to define the geometry, and its resolution. 
Many routines have been developed to generate grids within 
and about realistic surfaces that bound fluid flows. The 
computer program used to develop the grids used in this 
investigation. Turbomachinery C-Grid (TCGRID) , generates 
three-dimensional C- or K-type volume grids for 
turbomachinery. Grid files are compatible with RVC3D and 
PLOT3D [Ref. 7] which were used extensively in the flow 
solution and visualization. 

TCGRID generates a three-dimensional computational 
volume grid in the following manner: 

1. A course, equally - spaced meridional grid is generated 
between the supplied hub and tip. 

2 . The blade geometry is interpolated onto the meridional 
grid . 

3. 2-D blade - to -blade grids are generated along the 
meridional grid lines in (m, rbar* theta) coordinates 
using an old version of the Sorenson GRAPE code 

[Ref . 8] . 

4. The (m, rbar * theta) coordinates are transformed to 
(z,r, theta) coordinates. 

5. The 2-D grids are reclustered spanwise to make a full 
3-D grid. 
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6. Finally, the (z,r, theta) coordinates are transformed 
back to (x,y,z) and written in PL0T3D format. 

Hub, tip and blade geometries must be in MERIDL [Ref. 9] 
format . 

The hot -flow blade geometry data were received from NASA 
Lewis Research Center [Ref. 10] in (z,r, theta) coordinates 
appropriate for use in TCGRID . Since the Naval Postgraduate 
School facility was cold-flow, a 99% thermal - shr ink factor 
had to be applied to all linear dimensions. The FORTRAN 
code written to accomplish this may be seen in Appendix D. 

Grid generation and sizing is a subjective exercise in 
which there is no single ’right' answer. Often the success 
of a grid remains unknown until it is used with a flow 
solver . In this investigation, over one hundred flow 
solutions were attempted with various grid sizes before a 
successful combination was found. It was understood that 
the investigation required a fully three-dimensional viscous 
flow solution. Because of the resolution needed in the 
boundary layer, the grid had to be very fine. However ,i the 
successful grid contained more than twice the number of grid 
points in the spanwise direction than originally thought. 

The final grid size used for the flow solutions was 151 x 31 
x 45. Chima [Ref. 11] recommended 15 grid points in each 
endwall boundary layer for good resolution. An example of 
the TCGRID input file may be found in Appendix D. 
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Figures 8 through 11 show C-grids developed by TCGRID for 
the hub, midspan and tip. 

This geometry was difficult to model because of the 
endwall contour (Fig. 8) . The hub and tip geometry was 
defined in tne input file with (z,r) coordinates. Then the 
stator contour was defined in five sets of (z, theta, r) 
coordinates . Each set of coordinates was at a constant 
radius. Since the endwall radius reduced at the stator tip, 
some of the coordinates were defined outside the endwall. 
TCGRID calculates the intersection of the stator coordinates 
with the endwall and generates the grid defined by the 
intersection contour. Figures 9 through 11 show grids 
developed for hub, mean and tip sections. Figure 12 shows a 
three-dimensional view of the final grid. The curved 
trailing edge shown in this figure also complicated the grid 
generation procedure. 

B. FLOW SOLUTION 

1 . RVC3D 

Once a proposed grid has been developed, it can be 
used by a flow solver to simulate a flow. RVC3D (Rotor 
Viscous Code 3-D) is a computer code developed at NASA Lewis 
Research Center for the analysis of fully three-dimensional 
viscous flows in turbomachinery. The code solves the thin- 
layer Navier • Stokes equations with an explicit finite- 
difference technique. It is applicable to annular blade 
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Figure 10 . C-Grid at Stator Midspan 
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Figure 12. Three-Dimensional Stator Grid 
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rows or linear cascades and includes two turbulence models 



along with a simple tip clearance model. Mathematical 
formulation of RVC3D is described in references 12 and 13. 

The code solves the Navier - Stokes equations in 
Cartesian coordinates with rotation about the x-axis. The 
equations are mapped to a general body-fitted coordinate 
system (Fig. 13) . Turbulence is modelled using a three- 
dimensional adaptation of the Baldwin- Lomax model [Ref. 14] 
or the Cebeci -Smith model [Ref. 15] . A sample RVC3D input 
file may be seen in Appendix D. 

2 . Solution Method 

Simulations were run on two different computer 
systems; a Silicon Graphics Iris workstation and a Cray Y-MP 
EL-98 supercomputer. Both systems reside at the Naval 
Postgraduate School. Solutions were obtained for two inlet 
velocity profiles on each machine. The same version of 
RVC3D was loaded onto the two computers and identical input 
files were used. 

First, using the two computers for time efficiency, 
an inlet velocity profile to simulate a fully developed 
turbulent flow was found using the Iris while, concurrently, 
a profile simulating a slug flow inviscid profile was found 
using the Cray. The profiles were developed by modifying 
the RVC3D input files and running 100 iterations at a time 
to find the inlet velocity profile. After each run, the 
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Figure 13. General Body-Fitted Coordinate System 



31 



input file was modified and run again. The process was 
repeated until a satisfactory profile was obtained. It was 
assumed that the two computers would yield identical 
results . 

After several runs on each computer, it was noted 
that when similar values were used in the input files for 
the two machines, the velocity profiles that were calculated 
were vastly different. The turbulent velocity profile was 
run on the Iris until the solution converged and the same 
was done on the Cray with the slug profile. Then the input 
files were exchanged and the simulations were run again. 

All four solutions were plotted for comparison. 
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V. 



RESULTS AND DISCUSSION 



A . EXPERIMENT 

The data shown in Appendix C yield curves shown in 
Figures 14 and 15. The curves confirm general trends that 
were expected from the experiment. The data were taken 
without the intention of producing a turbine map but to 
validate the data acquisition methodology and to check the 
integrity of the turbine and its subsystems. The data show 
that the system will work but that it is inconsistent. It 
was noted that the water manometers used to measure pressure 
drop across the metering orifice were contaminated. This 
explains the inconsistent values noted in Appendix C as P31, 
P32 , P33 and P34 . 

The difference in the stagnation temperature readings 
indicate that the flow may not have been circumferentially 
uniform. The stagnation pressure values PT2 are obviously 
erroneous and the data obtained with that combination probe 
were not used. 

The low efficiencies shown in Figure 15 were expected. 
Since the shroud inserts were removed for these runs, it was 
expected that much of the air flow would pass between the 
blade tips and the outer case. The air passing around the 



33 



(s/uuqi) 0JBJJ MO|J SSB^ 




(uuqi/dH) JeMOd oypeds 



Figure 14. Specific Power and Mass Flow vs Speed 
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rotor blades did not contribute to the power output of the 
turbine and significantly reduced the efficiency. 

Figure 16 shows the speed sensor response of the TTR at 
800 rpm. The number in the lower left hand corner, measured 
in Hertz, comes from the rotational speed. But, since the 
magnetic speed sensor on the dynamometer detects the passage 
of teeth on a 30 tooth gear, the rpm of the system is twice 
the frequency shown. Figure 17 shows the speed sensor 
response at 7925 rpm. 

B. COMPUTATIONAL FLUID DYNAMICS 

As noted previously, running RVC3D with the same input 
files on the two computers resulted in solutions with 
clearly different inlet velocity profiles. Figure 18 shows 
the two turbulent inlet velocity profiles from one input 
file, and Figure 19 shows the slug flow profiles from the 
other. The differences in inlet profiles generated in the 
course of the solution, show that the inlet velocity profile 
is sensitive to either the architecture or the precision of 
the particular computer. 

Figures 20 and 21 show the differences in inlet velocity 
profiles using the two input files on the same computer . 

The similarity of the graphs shows that the output inlet 
velocity profile is not very sensitive to changes in input 
velocity profile. Comparison of Figures 18 and 19 with 
Figures 20 and 21 indicates that the output inlet velocity 
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Figure 16. Low Speed Vibration Response 
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Figure 17. High Speed Vibration Response 
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Figure 18 . Turbulent Inlet Velocity Profiles 
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Slug Inlet Velocity Profiles 
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Figure 19. Slug Flow Inlet Velocity Profiles 
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Normalized Velocity 



Cray Inlet Velocity Profiles 
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Figure 20. Cray Inlet Velocity Profiles 
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Normalized Velocity 



Iris Inlet Velocity Profiles 
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Figure 21. Iris Inlet Velocity Profiles 
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Normalized Velocity 



profile is more sensitive to computer than to input velocity 
profile . 

Figures 22 and 23 show the midspan chordwise variation 
in normalized pressure distribution of solutions with the 
same input velocity profiles on different computers. 

Figures 24 and 25 show the same results with the different 
inlet profiles on the same machine. Comparison of these 
four graphs shows that all the variations gave identical 
pressure distributions. Therefore, the pressure 
distribution is sensitive to neither inlet velocity profile 
nor computer . 

Figures 26 through 29 are PL0T3D plots of midspan Mach 
number contours. The first two show the Cray solutions of 
the turbulent and slug flow inlet velocity profiles. The 
latter two show Iris workstation solutions for the same 
input files. Again, in spite of the differences in the 
input files and the velocity profiles, all four solutions 
are identical. The plots show the flow accelerating on both 
the suction and pressure surfaces of the stator with maximum 
Mach number in the area of the throat as expected. Boundary 
layer growth as well as wake shape and orientation are 
evident in the plots. 

It is important to note that the wake, in yellow, does 
not follow the contour of the grid but follows the exit 
angle of the stator. The cut that joins the pressure and 
suction sides of the grid may be seen in Figure 10. The 
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Figure 22. Normalized Pressure - Turbulent 
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Figure 23. Normalized Pressure - Siug Flow 
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Figure 25. Normalized Pressure - Iris 
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Figure 26. Mach Contours - Cray - Turbulent 
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Figure 27. Mach Contours - Cray - Slug Flow 
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SSME HPFTP Turbine Stator Midspan Mach Contours 
Iris Solution - Turbulent Inlet Profile 




Figure 28 . Mach Contours - Iris - Turbulent 
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Figure 29. Mach Contours - Iris - Slug Flow 
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fact that the wake does not follow this cut indicates that 



the solution is somewhat independent of grid geometry and 
that continu'.ty is satisfied at each point in the flow. 

Figures 30 through 33 the total pressure contours across 
the exit pla re of the stator. These plots are also nearly 
identical bui. small variations do exist in the upper left 
corner of the exit plane. The turbulent solutions, Figures 
30 and 32, are different in this region from the slug flow 
solutions, Figures 31 and 33. This would indicate that the 
solution is somewhat dependent on inlet velocity profile. 

In these figures, it is also possible to see the 
variation in wake location along the span. The yellow area 
at the edges of the plots shows that the wake bows toward 
the pressure side at midspan. This may be explained by the 
bowed shape of the trailing edge of the stator . As the flow 
turns around the stator, the extended arclength along the 
surface in the region of the midspan may turn the flow 
farther than the shorter arclength along the stator at the 
hub and tip. 

In the boundary layer at the hub there is a small region 
of blue and yellow in the area of the grid cut. Although 
not a discontinuity, this locus coincident with the grid cut 
may indicate some solution dependence on grid geometry. 

An example of the residual histories of two solutions is 
shown in Figure 34 . The plot shows that the residuals 
matched exactly for the first 3500 iterations. The Iris 
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Figure 30. Total Pressure - Cray - Turbulent 
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Figure 31. Total Pressure - Cray - Slug Flow 
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Figure 32. Total Pressure - Iris - Turbulent 
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Figure 33. Total Pressure - Iris - Slug Flow 
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1000 2000 3000 4000 5000 6000 7000 8000 

Iteration 



solution was continued beyond 3500 iterations in order to 
see the effect of varying the inlet velocity profile. Each 
of the variations of the input file is marked by a 
significant rise in residuals. However, after approximately 
500 iterations each solution converged to nearly the same 
residual value. Even after seven changes in the input file 
to vary the inlet velocity profile, the residuals converged 
to four orders of magnitude below their start. This shows 
that the solution with the input grid generated by TCGRID is 
very stable. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



A. EXPERIMENT 

A working test model of the first stage of the SSME 
HPFTP ATD turbine has been put into operation in the Naval 
Postgraduate School's Turbine Test Rig (TTR) . The TTR 
subsystems have been successfully exercised; however, all 
lines and connections should be inspected, cleaned and 
repaired as necessary before any measurements are taken for 
performance evaluation. Connections should be completed to 
the water temperature probes and torque meter on the 
dynamometer and data acquisition should be automated. Three 
methods of calculating the power out of the turbine should 
be automated and output on-line: namely, using the flow rate 
and stagnation temperature drop of the air through the 
turbine, using flow rate and temperature rise of the water 
in the dynamometer, and using the rotational speed and 
torque generated in the dynamometer . Once the system 
alignment is verified, the shroud inserts should be 
installed and the performance measured. 

The inlet flow circumferential uniformity should be 
checked. If the flow is not uniform, a flow straightener 
should be designed and inserted into the inlet duct. It is 
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also recommended that the inlet and exit velocity profiles 
be determined experimentally . 

Ultimately, the turbine outer casing will need to be 
modified to allow LDV measurements of the flow, particularly 
in the tip clearance region. 

B. NUMERICAL SIMULATION 

A working CFD model of the SSME HPFTP ATD turbine first 
stage stator has been generated. Four solutions were 
generated on two different computers. It is recommended 
that more work be done to investigate the inlet -velocity 
profile dependence on the computer system, as this is not 
acceptable. An investigation is needed to determine whether 
the anomaly at the grid cut on the hub is grid dependent, or 
is a real, physical phenomenon. 

In order to yield more realistic flow solutions, the 
inlet struts that hold the bullet nose upstream of the 
first-stage stator should be modeled so that their effect 
can be included in the flow solution. Computations can be 
made with input of the actual inlet velocity profile, once 
this has been determined experimentally. For an analysis of 
the complete stage, the rotor geometry must yet be modeled. 
Such a modeling of the downstream portions of the turbine 
will allow comparison with experimental results when the 
turbine is tested with the rotor clearance at the design 
value . 
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APPENDIX A. ENGINEERING DRAWINGS* 
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Figure A1 . Bearing Outer Spacer 
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Figure A2 . Probe Mount 
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APPENDIX B . CALCULATIONS 



A. MASS FLOW RATE 

The following assumes the 9.0 in. orifice is installed 
and is derived from Vavra [Ref. 16] . 

First solve: 



Now find; 



where , 



w '= ca Y 1 C 



ifi'=ca Y 1 



N T i 
N T i 



with C=l; 



w 

x=0 .2525 — — 
z 



Z = 1.9+0.24 ( i--l) 

100 



C = f(x) 

in lbm/s. 

tj= temperature ahead of orifice ( C F) 

\.= pressure drop across orifice (in. H 2 0) 

P x = absolute pressure ahead of orifice (in. Hg) 
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1 . Flange Taps 



m'=ca Y 1 



M T 1 



where T 1 = t 1 + 460 
Therefore , 






rn'-S . 6445 [1*0 . (1-0 . 027 ±') 



z= 1.9+0.24 ( — — -1 ) 
100 



rb I 

x=0 . 2525 — 
z 

r_ lJL o . 0047 
x 



ri)=£rii / ’. 



2 . Vena Contracta Taps 



lh'=6 .7720 [1 + 0 . 0015 



(ti- 60 ) 

100 



] (1-0.027-!?) 

i i \ 



p iK 

^ + 460 



Z=1 .9+0.24 ( 



100 



- 1 ) 



if}^ 

x=0 . 2525 — 
z 

0 . 0029 

* x 



/n=(/n / . 
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B . TURBINE POWER 



Power =rfiC ( AT 

C p = 0.24 Btu/lbm.'R for air. 

AT=TT3-TT4 

TT4 = Stagnation Temperature of Air at Turbine Exit ( * F ) 

TT3 = Stagnation Temperature of Air at Turbine Entrance (’F) 
Since 1 Hp°s = 550 ft.lbf and 778 ft.lbf = 1 BTU ( 

Power=m • 0 . 24 ( TT3-TT4) 778 . 

550 



C. TURBINE EFFICIENCY 



Efficiency= 



TT 3 -TT 4 

TT 3 -TT 4i 



where TT 4i = Isentropic stagnation temperature is given by 

pt. n 

TT,, =TT, ( i) t . 

4 1 3 ' prji t 
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APPENDIX C 



COMPLETE TABLE OF MEASUREMENTS 
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APPENDIX D. FORTRAN PROGRAM AND INPUT FILES 



RVC3D Input File 



' SSME HPFTP -- Indigo -- Real Slug Flow' 

&nll im=151 jm=31 km=45 itl=20 iil=70 &end 

&nl2 cf 1=4 . 0 aviscl=0.0 avisc2=0.5 avisc4=1.5 ivdt=l nstg=4 
itmax = 9 85 

irs=l epi=0.50 epj=0.60 epk=0.60 &end 

&nl3 ibcin=3 ibcex=3 isymt=0 ires=l icrnt=50 

iresti=l iresto=l ibcpw=0 iqin=0 &end 

&nl4 emxx=0.212 emty=0.0 emrz=0.0 expt=0.0 prat=0.6500 
ga = l . 4 

om=0 . 000000 igeom=l alex=-60.0 &end 
&nl5 ilt=2 tw=l . 00 renr=6.651e6 pxnr=.7 prtx=.9 
vispwx= .666666 

sxtip=0.0 cmutrri=14 . jedge=15 kedge=ll iltin=0 dblh=0.00 
dblt=0 . 00 &end 

&nl6 iol=l io2=151 oai=0. ixjb=0 njo=l nko=3 
jo=l ko= 13 17 21 &end 
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Inlet. f - Reads RVC3D output and builds files for plotting 
Inlet Velocity Profiles. 

(2k k k k k k k k k k k * k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k 

c inlet . f reads rvc3d files & writes ascii files for 

plotxy 

c unit 1 = input xyz file (fort.l) 

c unit 3 = input q file (fort. 3) 

c unit 17 = output inlet velocity profiles 

c unit 4 = output residual history 

c*****-*****************************************-************* 

parameter (ni=152,nj=32, nk=46 ) 

real x (ni , nj , nk) , y (ni , nj , nk) , z (ni , nj , nk) 

real qq (5 , ni , nj , nk) , resd (5000 , 5) 

real vel (5) , zz (5) 

dimension ii(5) 

c i- values are hard-wired below 

data ii/60 , 70 , 75, 80, 90/ 

£j***********A*** ********************************************* 

c read grid coordinates 

Q***********-K *********************************************** 

r ead ( 1 , * ) im , j m , km 

read (1 , * ) ( ( (x (i , j , k) , i=l , im) , j =1 , jm) , k=l , km) , 

1 ( ( (y (i , j , k) , i = l , im) , j =1, jm) , k=l , km) , 

2 ( ( (z (i , j , k) , i = l , im) , j =1, jm) , k=l ,km) 

£*kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

c read restart file 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

r ead ( 3 , * ) imax , j max , kmax 
read (3 , * ) f smach, alpha , re, time 
c 

icheck=iabs (im- imax) +iabs ( jm- jmax) +iabs (km- kmax) 
if ( icheck . ne . 0 ) then 

write (6 , 6 10 ) im, jm, km, imax, jmax , kmax 
stop 
endif 
c 

read (3 ,*)((( (qq (1 , i , j , k) , i = l , im) , j =1 , jm) , k=l , km) ,1 = 1,5) 
c 

c additional residual data 

read (3,*)itl,iil, phdeg , ga , om, nres , dum, dum, dum, dum 
read (3 , * ) ( (resd (nr , 1) , nr=l , nres) ,1=1,5) 

Q kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

c vel to unit 17 

Qk k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k 

itr=im+l - itl 
ggm=ga* (ga- 1 . ) 
j= jm 

c normalize x by chord 



68 



c do 7 1-1,5 

c k=kk ( 1 ) 

c xmin=x ( im/2 , j , k) 

c xmax=x(itl, j,k) 

c do 5 i=itl,itr 

c xmin=min (xmin, x ( i , j , k) ) 

c 5 xmax=max (xmax , x ( i , j , k) ) 
c chord=xmax -xmin 

c do 7 i=itl , itr 

c 7 x ( i , j , k) = (x ( i , j , k) -xmin) /chord 
c 

i=7 5 
j=jm 
k=km/2 

velinf = sqr t ( (qq ( 2 , i , j , k) * *2+qq (3,i,j,k)**2 
1+qq (4 , i, j , k) **2) /qq ( 1 , i , j , k) ) 
c 

c write(7,305)(kk(l), 1=1,5) 

do 20 k=l,km 
do 10 1=1,5 

i=ii (1) 

vel ( 1 ) =sqr t ( (qq (2 , i , j , k) * *2+qq (3 , i , j , k) * * 2 
1+qq (4 , i , j , k) * * 2 ) /qq (1, i, j , k) ) 

10 zz(l)=z(i,j,k) 

20 write(17,300)i, (zz(l) , vel (1) /velinf, 1=1, 5) 

£-*•*■*•*★★■*■*■*•*•*■*■*•*•*•*■*■*■***■**•**•*•*•*■*•*■*■*•*■*:**•*■*■*■*•*•*■*•*■*•**•*•** 'k'k’k'k'k'k'k + ic 

300 format (i5 , lOf 10 . 5) 

305 format ( ' k= ' , 5 ( 17x , i3 ) ) 

610 format ( ' ***** warning *****',/, 



1 

match ' , / , 


' im, 


jm, 


km. 


read 


from 


input', 3i5,’ do not 


2 

stop 

end 


' im, 


jm, 


km, 


read 


from 


restar t f ile ' , 3i5) 
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program shrink.99 

This inputs data and multiplies it by 0.99 for a 99% 
shrink factor . 



real x (988) ,p (988) 

open (uni t = 2 , f ile= ' hotlv.dat' ,status='old' ) 
open (unit=15, file= ' lv. tcgin' , status= ' new' ) 
read(2, A ) (x(i),i = l,988) 
do 20 i =1 , 9 88 

p (i) =x ( i ) * 0 .99 
20 continue 



print 


the 


data out in blocks 


write 


(15, 


* ) 


(p (i) , i=l , 88) 


wr i te 


(15, 


*) 




wr i te 


(15, 


M 


(p ( i ) , i=89, 176) 


wr i te 


(15, 


M 




wr i te 


(15, 


* ) 


( p ( i ) , i=177 , 264) 


wr i te 


v 15 , 


* ) 




wr i te 


(15, 


* ) 


(p (i) , i = 26 5 , 352) 


wr i te 


(15, 


* ) 


■ 


wr i te 


(15, 


A ) 


( p ( i ) , i=353, 405) 


write 


(15, 


A ) 




wr i te 


(15, 


* ) 


(p ( i) , i = 4 06 , 4 58) 


wr i te 


(15, 


A ) 




wr i te 


(15, 


A ) 


( p ( i ) , i=459 ,511) 


wr i te 


(15, 


A ) 




wr i te 


(15, 


A ) 


(p ( i ) , i=512, 564 ) 


wr i te 


(15, 


A ) 




write 


(15, 


A ) 


(p ( i ) , i = 56 5 , 6 17 ) 


wr i te 


(15, 


A ) 




wr i te 


(15, 


A ) 


(p ( i ) , i = 6 18 , 67 0 ) 


wr .i te 


(15, 


A ) 




wr i te 


(15, 


A ) 


(p ( i ) , i =67 1,723) 


wr i te 


(15, 


A ) 




wr i te 


(15, 


A ) 


(P ( i ) , i=724 ,776) 


wr i te 


(15, 


A ) 




wr i te 


(15, 


M 


(p ( i ) , i=777 , 829) 


wr i te 


(15, 


A ) 




wr i te 


(15, 


A ) 


(p ( i ) , i=830, 882) 


wr i te 


(15, 


A ) 




wr i te 


(15, 


A ) 


(p ( i ) , i=883 , 935 ) 


wr i te 


(15, 


A ) 




write 

end 


\ 15 , 


A ) 


( p ( i ) , i =9 36 , 9 88 ) 
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TCGRiD Input File 





'i 

<D 

n. 


.^-rM.nnnnM.w,,. .., -T 7 ? V 'i ' VTi ** 




oSJXj:" ^*11 s ,'o» s.?,?; f’s^r;',.**; ^r,r? !: , sr;*t! t*-* ?:*■“••; -:■ ; •. ; ..,««. r rr ~*-„ . rrr ~ 


e 

—7 


r V . i * i r < 


s 


KaassasaassKaasasaaaar.fiS^iasssassr::.;;;;;; K:.:-;:.."-..3:ss-2SSSs:jsss?r.j 

a si 3 a a a ;:; 3 s a 3 3 3 a a a 3 a a a 33n3Ssr»«£;; a s r 5 « ■! r. H'-.iiSftSniVi « i H ?• ?; H 3 3: • • • •" -; :-, rFrrr.r .-■ F P R » 


■ ■ ‘ :f.pr?'r ; ' ' • --■•-• ■■•■ ' : •- •••' ; 

aaaaaaaaaaaiiaaaHa^ 

aaKa3333a3338a3a3aaa?JiH’t«Sr35if533:'s« : ir.3 5 ,3:;3n3a:irjRa?.a333h?/i}n:;nr ! *???rrrrPP 


U) 

g 




o 

ro 


(NM^MMNNm<niir(NfMryiMfMNrNHi»('(B»',i,ir<nv^iMr v» -» ■ ■ < hi il j <i* rr 11 «v *-i < i < > oi . . i/i i . 4 r i •*■ i in ^ r. . . t r- <n «■ <tt cr aictn®oia)<i)®'r 

inui.ri.iiiiKiin iniruiiiiuiiih/uoinim/U'Hii » . n w mo^ n m»>Ki »n\f h 'P i • \» n . i m w ’o vf n » m in n >n o m r- O if m a irffnin'P'crwa'a 

n sss ^ srjrj:*.r.^r 


& 

1 


- ?777?7 '• .Vi 




7 2 





7 3 
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